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wedge  is  modeled  as  a  point  reacting  complex  quantity  which  is  independent  of 
the  other  face.  The  solution  was  constructed  as  an  angular  spectrum  to  satis¬ 
fy  the  boundary  conditions  and  Sommerfeld  radiation  condition.  The  solution 
kernel  was  obtained  exactly  and  is  in  terms  of  circular  functions.  The  sol¬ 
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by  use  of  asymptotic  techniques.  This  solution  is  much  simpler  than  the  one 


00  ,  :?nRM7,  1473 


EOITION  OF  I  NOV  6S  1$  OBSOLETE 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  RAGE  flWian  Oat*  Enttt+r' 


Mi, 


( 


MCUWTy  CUAMiriCATIOM  OF  THIS  »AQ«flW>i  Dim  KmUml) 

developed  by  Russian  scientists  (  for  example,  see  G.D.  Maliuzhinets ,  "The 
Radiation  of  Sound  by  Vibrating  Boundaries  of  an  Arbitrary  Wedge,"  Farts 
1  and  2,  Soviet  Physics  Acoustics,  pp.  152-174  and  240-248  (1955)  )  which 
was  obtained  by  a  method  similar  to  Wiener-Hopf  technique.  Thus,  it  is 
easier  to  use  in  highway  noise  applications  because  of  its  simplicity. 

The  solution  for  the  diffracted  pressure  exhibits  clearly  the  role  of 
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The  prediction  of  the  sound  scattered  by  impedance  covered  wedges 
is  obtained  by  use  of  dual  integral  equations.  The  impedance  of  each 
face  of  the  wedge  is  modeled  as  a  point  reacting  complex  quantity  which 
is  Independent  of  the  other  face.  The  solution  was  constructed  as  an 
angular  spectrum  to  satisfy  the  boundary  conditions  and  Sommerfeld 
radiation  condition.  The  solution  kernel  was  obtained  exactly  and  is 
in  terms  of  circular  functions.  The  solution  of  the  scattered  pressure 
was  then  obtained  for  far-field  and  mid-range  by  use  of  asymptotic 
techniques.  This  solution  is  much  simpler  than  the  one  developed  by 
Russian  scientists  [for  example,  see  G.D.  Maliuzhinets ,  "The  Radiation 
of  Sound  by  Vibrating  Boundaries  of  an  Arbitrary  Wedge,"  Parts  1  and  2, 
Soviet  Physics  Acoustics,  pp.  152-174  and  240-248  (1955)1  which  was 
obtained  by  a  method  similar  to  WLener-Hopf  technique.  Thus,  it  is 
easier  to  use  in  highway  noise  applications  because  of  its  simplicity. 

The  solution  for  the  diffracted  pressure  exhibits  clearly  the  role 
of  the  incident  and  reflected  shadow  boundaries  and  shows  there  is  one 
minimum  in  the  scattered  field  which  depends  on  the  two  surface 
impedances.  For  backscattered  pressure,  the  solution  exhibits  two 
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CHAPTER  I 


1.2  Literature  Review 

The  problem  of  diffraction  by  sharp-edged  objects  has  been  studied 
extensively  in  acoustics,  electromagnetics,  optics  and  water-wave 
theory.  The  following  literature  review  is  mainly  concentrated  on 
diffraction  by  wedges. 

The  first  investigator  who  succeeded  in  deriving  the  first  exact 
solution  of  diffraction  of  a  plane  wave  by  a  wedge  was  Sommerfeld  [1]. 
His  solution  is  only  applicable  to  a  particular  case  of  two-dimensional 
problem  of  diffraction,  namely  two  types  of  ideal  boundary  conditions 
where  both  surfaces  are  either  perfectly  soft  (pressure  release)  sur- 
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faces  or  perfectly  hard  (rigid)  surfaces.  The  Sommerfeld  theory  has 
been  considered  to  be  very  difficult  to  understand  because  the  solu¬ 
tion  is  derived  by  a  heuristic  image  method.  Sommerfeld  arrived  at  his 
results  by  seeking  an  appropriate  solution  to  the  wave  equation  of 
period  4ir  and  by  combining  it  with  its  "image."  The  solution  that 
resulted  is  valid  in  the  far-field,  but  it  becomes  unbounded  at  the 
incident  and  reflected  shadow  boundaries.  MacDonald  [2]  approached  the 
wedge  diffraction  problem  using  the  classical  separation  of  variables 
technique  and  expressed  the  solution  as  an  infinite  series  of 
appropriate  eigenfunctions  for  point  and  line  sources. 

The  first  solution  of  diffraction  by  an  imperfectly  conducting 
obstacle  was  obtained  by  Senior  (3].  The  general  problem  of 
diffraction  of  a  plane  wave  by  a  non-ideal  surface  wedge  of  arbitrary 
angle  was  considered.  In  order  to  obtain  the  solution.  Senior  made 
some  assumptions  without  physical  reasoning  which  hold  only  for  a 
pressure  release  wedge.  Williams  [4]  obtained  a  solution  for 
diffraction  by  impedance  covered  arbitrary  wedge.  He  reduced  the 
original  problem  to  the  solution  of  an  ordinary  difference  equation. 
This  equation  was  then  solved  in  terms  of  the  double  Gamma  function 
defined  by  Barnes  (5j.  The  resulting  solution  is  an  infinite  product 
form  and  can  only  be  solved  in  closed  form  if  the  wedge  angle  is  np/q 
where  p  and  q  are  relative  prime  integers,  with  p  odd. 

Russian  physicist  Maliuzhinets  [6]  demonstrated  chat  the  kernel 
of  the  problem  obtained  in  the  Fourier  transform  approach  is  the 
solution  of  a  difference  equation.  He  gave  conditions  and  proved 
theorems  about  this  fact.  However,  his  method  for  obtaining  the 


solution  of  this  difference  equation  was  not  published  in  the  open 
literature.  The  solution  is  quite  similar  to  the  Weiner-Hopf  tech¬ 
nique  of  factoring  a  function  into  the  product  of  two  functions  each 
valid  in  one  half  space.  Skudrzyk  [7]  and  VanLennep  [8]  made  two 
independent  attempts  and  derived  the  so  called  Maliuzhinets  functions. 
Maliuzhinets  functions  consist  of  the  product  of  four  complicated 
functions  with  four  different  arguments.  These  functions  are  typical 
of  those  found  in  Weiner-Hopf  solution,  i.e.,  they  are  exponential 
functions  each  of  whose  argument  is  a  proper  integral  of  another 
function.  Several  applications  of  these  functions  to  propagation 
problems  were  published  by  Maliuzhinets  [9-11]  and  other  Russian 
scientists  [12-15].  The  closed  form  solution  obtained  by  using 
these  functions  is  limited  by  the  same  conditions  applied  to  the 
solution  of  Williams.  For  a  general  wedge  angle  the  solution  was 
approximated  by  Skudrzyk  [ 7 ] . 

In  more  recent  studies  Pierce  [16]  and  Kouyoumjian  [171  extended 
the  available  solutions  of  Williams  and  Sommerfeld  for  applications 
in  electromagnetic  waves  and  acoustic  absorbing  barriers,  respectively. 
Pierce  was  interested  in  the  application  to  more  complicated  barrier 
structures  which  involve  more  than  one  wedge  —  thick  absorbing 
barrier  or  trapezoid  absorbing  barrier  —  by  using  the  diffraction 
coefficient  for  the  absorbing  wedge  in  conjunction  with  Keller's 
Geometric  Theory  of  Diffraction  (GTD) .  In  Kouyounjian’s  original 
paper  only  diffraction  by  a  perfectly  conducting  wedge  whose  edge 
is  a  curved  line  was  studied.  The  diffraction  coefficient  obtained 
by  incorporating  asymptotic  results  by  Pauli  [18]  and  Oberhettinger 
[19]  is  continuous  across  the  incident  and  reflection  shadow  bound- 
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aries.  The  expressions  for  the  acoustic  wedge  diffraction  coeffi¬ 
cients  contain  Fresnel  integrals,  which  ensure  that  the  total  field 
is  continuous  at  the  incident  and  reflection  shadow  boundaries. 

These  diffraction  coefficients  were  modified  for  non-ideal  boundary 
conditions  by  using  Maliuzhinets  functions.  A  complete  summary  of 
these  solutions  were  published  by  Rayek,  et  al.  [20].  Rawlins  [21] 
proposed  a  superposition  method  for  perfectly  absorbing  wedges  by 
combining  the  known  solutions  for  hard  and  soft  wedges  and  derived 
a  simple  closed  form  solution  for  an  absorbing  wedge. 

The  most  widely  studied  special  wedge  is  the  half-plane. 

Carslaw  [22]  extended  Sommerfeld’s  plane  wave  solution  for  a  half¬ 
plane  with  ideal  (either  rigid  or  soft)  boundary  conditions  for  line 
and  point  sources.  MacDonald  [23]  showed  that  his  wedge  solution 
reduces  to  Carslaw’ s  half-plane  solution.  He  also  derived  asymptotic 
solutions  in  terms  of  Fresnel  integrals  which  are  widely  used  in 
modern  diffraction  analysis.  Rawlins  [24]  obtained  a  closed  form 
solution  for  diffraction  by  a  semi-infinite  plane  with  rigid-soft 
surfaces  by  using  a  modified  Wiener-Hopf  method  (which  is  a  standard 
method  for  half-plane  problems). 

The  first  exact  solution  of  diffraction  by  a  semi- infinite 
metallic  sheet  with  finite  conductivity  was  obtained  by  Senior  [25]. 
He  considered  this  particular  case  of  wedge  by  means  of  the  Wiener- 
Hopf  technique.  The  result  is  in  integral  form  and  cannot  be  easily 
computed.  Clemnow  [26]  introduced  a  "dual  integral"  method  for 
ideal  boundaries.  In  a  more  recent  studv  Kendig  and  Havek  [27] 
employed  an  extension  of  a  dual  integral  method  for  a  hard-soft 
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barrier.  Kendig  [28]  also  utilized  the  sai_e  method  to  solve  for  the 
diffraction  by  an  impedance  covered  half-plane.  Kendig  and  Hayek 
obtained  by  Fourier  transforms  two  integral  equations,  each  holding  for 
a  different  side  of  the  half-plane.  These  dual  integral  equations  were 
solved  and  continuous  solutions  at  the  incident  and  reflection  shadow 
boundaries  for  the  plane  wave,  line  and  point  sources  were  derived. 

Extensive  reviews  of  diffraction  from  wedges  are  given  by  Bowman, 
Senior,  Usleghi  [29],  Pierce  [16],  Christiansen  [30],  Skudrzyk  [7]  and 
Hayek,  et  al.  [20]. 


CHAPTER  II 


DIFFRACTION  OF  PLANE  WAVE  BY  AN  IMPEDANCE  COVERED  WEDGE 

2.1  Statement  of  the  Problem 

When  an  acoustic  wave  is  incident  upon  a  discontinuous  body, 
having  acoustic  characteristics  different  from  those  of  the  surrounding 
medium,  the  resulting  acoustic  field  is  modified  from  the  incident 
wave.  The  modified  wave  field  is  caused  by  scattering  from  the  con¬ 
tinuous  and  discontinuous  parts  of  the  body.  Therefore,  in  addition  to 
the  initial  field,  there  is  scattered  field  which  is  propagated  outward 
away  from  the  discontinuity  and  interfering  with  the  incident  field. 

In  this  study,  solution  is  sought  for  the  diffraction  of  an 
incident  plane,  line  or  point  source  by  a  wedge  with  an  arbitrary  wedge 
angle  and  covered  by  surfaces  possessing  locally  reacting  impedances. 

It  should  be  pointed  out  that  even  diffraction  from  a  half-plane,  which 
is  one  limiting  case  of  a  wedge,  is  approached  by  ideal  boundary 
conditions.  These  ideal  surface  conditions  are  characterized  as  either 
the  pressure  release  or  the  rigid  cases.  In  this  study  more  realistic 
and  mathematically  more  complicated  surface  conditions  are  to  be 
considered.  The  acoustic  impedances  of  the  wedge  surfaces  are  finite, 
locally  reacting  and  complex.  The  real  and  imaginary  parts  of  the 
complex  impedance  are  the  acoustic  resistance  and  reactance, 
respectively.  The  acoustic  resistance  of  the  surface  impedance 
represents  the  acoustic  energy  that  is  absorbed  by  the  surface  of  the 
wedge.  The  acoustic  reactance  of  the  surface  impedance  represents  the 
Inertia  and  the  compliance  of  the  surface  coverage. 
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One  can  conceive  a  surface  possessing  .ocally  reacting  point 
impedance  Z  on  which  p*Zvn,  with  vn  being  me  normal  velocity  and  p  is 
the  surface  pressure  holds  at  each  and  every  point.  This  impedance  is 
independent  of  the  angle  of  incidence  of  the  acoustic  wave.  The 
locally  reacting  point  impedance  model  works  well  for  materials  such  as 
grass  covered  ground,  porous  media,  fiberous  material  such  as 
fiberglass  wool,  etc.  However,  this  model  ignores  the  wave  effects  on 
the  surface  or  the  effects  of  refraction  in  the  material  of  the  wedge. 

A  locally  reacting  surface  can  be  best  visualized  as  a  surface  which 
consists  of  isolated  simple  oscillators.  In  other  words,  there  is  no 
coupling  with  the  neighboring  points  of  the  surface,  which  excludes 
elastic  waves  propagating  on  the  surface. 

2.2  Basic  Approach 

The  solution  of  the  problem  is  partially  based  on  the  works  of 
Clemmow  [26],  Kendig  [28],  and  Kendig  and  Hayek  [27].  It  begins  with 
writing  separately,  the  equations  representing  the  acoustic  field 
generated  by  a  plane  wave  source  valid  in  the  two  half  spaces  symmetric 
about  an  infinite  wedge.  By  employing  Fourier  transform  techniques  and 
taking  advantage  of  the  symmetry  of  the  wedge,  two  integral  equations 
are  obtained.  The  solution  was  constructed  as  an  angular  spectrum 
function  which  satisfies  these  two  integral  equations  and  the 
Sommerfeld  radiation  condition.  Since  the  acoustic  field  can  be 
produced  by  various  kinds  of  sources,  the  scattered  field  due  to  line 
and  point  sources  is  also  investigated.  The  scattered  pressure  was 
Chen  obtained  for  far-field  and  mid-range  by  use  of  asymptotic 
techniques. 


2.3  Integral  Representation 

The  wedge  is  defined  in  terms  of  cylindrical  coordinates 
(r,4>,z)  by  the  equations  4>*8  (upper  surface)  and  (lower 

surface)  and  its  edge  is  coincident  with  the  z-axis  as  shown  in 
Fig.  2.1. 

A  scattering  problem  consists  of  finding  the  scattering 
pressure  p  when  an  incident  plane  wave  p^  impinges  on  the  obstacle. 
The  total  field  pressure  must  satisfy  three  requirements,  the 
differential  equation,  the  boundary  conditions  and  a  continuity 
condition.  The  standard  Fourier  method  is  not  applicable  to  a 
space  with  a  discontinuity.  Fourier  integrals  converge  only  in 
half-space.  To  assure  the  continuity  at  the  intersection  of  two 
half  spaces,  the  continuity  of  the  pressure  and  the  normal  velocity 
are  imposed. 

Consider  the  homogeneous  wave  equation  for  pressure 


c2  3t2 


(2.1) 


where  c  is  the  sound  speed  of  the  acoustic  medium 

Assuming  the  pressure  to  be  periodic  in  time  as  p(x,t)  = 
p(x)  then  Eq.  (2.1)  reduces  to  the  Helmholtz  equation: 


'p  +  k'p  =  0 


(2.2) 


where  k  *  j/c 


10 


Define  the  Fourier  Transform  pair 

(2.3) 

(2.4) 

Applying  the  above  transform  on  x  in  Eq.  (2.2),  one  obtains: 


f(x)  •  i:  f  F(y)  e'lYXdy 

"  _ao 

CD 

F(y)  *  /  f (x)  eiyxdx 
-.00 


Since  the  wave  equation  is  of  the  second  order,  it  has  two 
independent  solutions.  One  solution  has  to  be  discarded  because 
it  represents  energy  sources  at  infinity.  The  diffracted  pressure 
is  directed  away  from  the  wedge  boundaries,  so  the  solution  must 
describe  outgoing  waves.  Using  plus  sign  (+)  for  y  <  0,  and  the 
minus  (-)  sign  for  y  >  0,  the  requirement  of  waves  propagating 
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and  decaying  away  from  the  wedge  is  satisfied. 

Using  the  inverse  transform,  Eq.  (2. J) ,  solution  of  the  wave 
equation  is  obtained  as 

t  \  k  (  .  ,  .  ±iky/l-A2  -ikAx,. 

p(x,y)  -  -zz  f  A(y)  e  3  e  dA 

— OO 


where  A  ■  Y/k,  and  A(X)  is  the  spectral  amplitude  function.  The 
above  integral  represents  the  solution  of  the  wave  equation  when 
integrated  in  the  proper  convergent  space. 

Although  the  convergence  properties  of  the  integrand  are  not 
affected  by  a  coordinate  transformation,  at  this  point  it  is 
advantageous  to  introduce  polar  coordinates. 

Let 


x  *  rcos0 
y  *  rsin<J> 


(2.9) 


and 


A  *  cosa 

The  exponent  in  the  integrand  thus  becomes 


-ikr(cosO  cosa±sind  sini)  (2.10) 

=  -ikrcos (£+a) 

The  limits  on  A(-°°  to  +»)  are  given  by 

A  =  cosj.  *  cos  (a  ) 

(2.11) 

■  cosct^  coshou  -  isin.xr  sinh 
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where  subscripts  r  and  i  designate  real  ard  imaginary  parts.  Due 
to  periodicity  of  cosa,  the  entire  A-spac-  is  mapped  into  a  strip 
of  7T  width  which  repeats  with  periodicity  Certain  values  for 
ar  and  must  be  chosen  in  order  to  determine  the  mapping  of  the 
A-path  in  the  orspace. 


ever ^  ?  su“'*s 


-r.  /■■  »■■  -0*  -C^-TV 
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Re[-ikr{cos(^-ar)coshai  +  isinha^  sir^-ar>}] 

(2.15) 

*  krsinha^  sin(4>-ar)  <  0 

If  <  0,  then 

<p  +  (2n-l)fT  <_  ar  <  2nTT  +  <p  (2.16) 

If  >  0,  then 

<)>  +  2nfr  ar  <_  <j>  +  (2n+l)7t  (2.17) 

Similarly  for  y  <  0: 

If  <  0,  then 


<j>  +  (2n-l)TT  +  2niT 


If  >  0,  then 


(2.18) 


-<p  +  2n7T  <  ar  <  -<p  +  (2n-l)T  (2.19) 

It  is  therefore  necessary  for  the  convergence  of  the  integral 
that  the  extremes  of  any  path  of  integration  must  lie  in  the  shaded 
regions  of  Figs.  2.3a  and  2.3b. 

The  part  of  the  path  along  the  real  axis  represents  homogeneous 
plane  waves  fanning  out  to  free  space,  whereas  the  vertical  parts 
represent  inhomogeneous  waves  decaying  away  from  the  wedge.  The 
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spectrum  function  P(cosa)  specifies  in  terms  of  amplitude  and  phase 
the  weight  attached  to  each  plane  wave  of-  the  spectrum.  By  this 
physical  reasoning  the  above  integral  could  have  been  written  down 
heuristically  as  superposition  of  plane  waves  propagating  in  all 
directions  multiplied  by  a  weighting  function  which  contains  the 
amplitude  and  the  phase  of  each  plane  wave. 

To  characterize  the  scattered  wave  caused  by  the  wedge,  it 
must  be  required  that  the  scattered  field  does  not  have  a  source 
at  infinity.  This  condition  is  described  by  the  "Sommerfeld 
radiation  condition." 

lim  {*7  +  ikp]}  -*•  0  (2.20) 

r  -*•  oo 

Substituting  the  expression  for  scattered  pressure  into  the 
above  expression  yields 

lim  (.7  [/  P(cosct)(l-cos(^))e'ikrcos(:+l)d-t]}  -  0 

r  -  oc  (2.21) 

This  limit  is  satisfied  if  the  path  of  integration  is  in 
the  convergence  regions  of  the  i-plane,  which  are  shown  as  the 
shaded  regions  in  Figs.  2.3a  and  2.3b.  So  the  right  choice  of 
y-paths  that  fall  in  the  convergence  regions  will  automatically 
satisfy  the  Sommerfeld  radiation  condition. 
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2.4  Plane  Wave  Formulation 

An  incident  plane  wave  in  free  space  can  be  represented  as 
follows : 


p.  =  p  e 
*1  *o 


ikrcos  ($-<{>  ) 
o 


(2.22) 


where  pQ  is  the  amplitude  of  the  incoming  wave  and  <pQ  is  the  angle 
the  wave  normal  makes  with  the  positive  x-axis,  see  Fig.  2.1. 

The  surface  impedance  is  the  ratio  of  the  total  pressure  to 
the  normal  component  of  the  velocity.  Then  the  boundary  conditions 
that  have  to  be  satisfied  are 


(2.23) 


where  pfc  is  the  total  pressure  and  vtn  is  the  normal  component  of 
the  total  field  particle  velocity. 

The  reflection  coefficient  for  a  plane  wave  incident  at  a 
plane  surface  with  impedance  Z  is  given  by 


cosu  -  -~c/Z 
cosu  +  Cc/Z 


(2.24) 


where  u  is  the  angle  of  incidence  with  the  normal  to  the  surface. 
There  will  be  no  reflection  if  the  numerator  is  zero,  that  is  the 
angle  of  incidence  is  equal  to  the  grazing  acoustic  Brewster  angle 


given  by 


_ ......  ■  i'>  ■ 

* r* w .% .  ■  -*  •  *  -•  ■  •  * 


v  -;•  ,k  %“V  v  V 


sin9+  =  ~ 
Z 


(2.25) 


(2.26) 


Then  the  boundary  conditions,  Eq.  (2.23),  can  be  written  as 


.+  _  l/°ptotal. 


ikp_  .  ,sin0  *  — (-srr 
“total  r  aip 


(2.27) 


ikptotalsin8' 


l.8p total. 
"  rh<p  > 


d>=  -B 


(2.28) 


or  in  more  familiar  form 


total 


1  3p - +  ikrsitl0~p total  =  0  iS 


(2.29) 


Using  the  expressions  found  in  Eqs.  (2.22)  and  (2.14)  for 
the  incident  and  the  scattered  pressure  the  boundary  conditions 
Eqs.  (2.29)  for  5*3  and  $  »  -B,  respectively  become 


/,.P(  cos:.)  [sine*  +  sin(3-:.)  Je-ikrcos(3_l)d. 


-?  (sin.+  -  sin(S-£  ))eikrCos(£":o) 


(2.30) 
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/  P(cosoi)  [sin9  +  sin(0+ot)]e  i^rcc- (8+ot)^ 


■  -p  (sino  -  sin(3+6  )) J-^-rcos 
o  c 


(2.31) 


For  reasons  which  will  be  explained  in  the  next  section  it  is 
beneficial  to  modify  the  boundary  conditions  into  a  different  form. 
First  consider  the  left  hand  side  of  the  boundary  condition  for 
<p  m  S,  Eq.  (2.30),  on  two  different  Y-paths,  Y^  and  Yj,  shown  on 
Fig.  2.4. 

/7T+i  P(cosa) (sin3+  +  sin(S-ct) ]e-i^rcos '^d  < 

-i°° 

ikrcosj 

-  /  P(S-a-ir)  [sin3  -sina]e1KrcoSadi  (2.32) 

-2~+6-i°° 


and 


j~2~~  i°° 
-T+i» 


P(cosa) [sino  + 


sin(3-‘.t)  ]e 


-ikrcos(3-  0, 


,£-i»  _»  r  .  .+  ...  ikrcos.i. 

/  P(3*>')[sui"  -  single  d.i 

-+2+^ 


(2.33) 


where  a  new  variable  :«.*  *  2-  --  is  introduced  in  both  of  the  above 

integrals,  i*  is  replaced  by  t  and  the  limits  of  integration  are 
changed  correspondingly.  Adding  the  expressions  in  Eqs.  (2.32) 
and  (2.33),  a  new  path  is  established  as  shown  in  Fig.  2.5,  so 
that  Eq.  (2.30)  reduces  to 


,  v  ,  .  .+  .  ,  ikrcosa. 

f  P(B-ct-iT)  [sin0  -  sina]e  da 

r 

-  -2p  (sin9+  -  sin(3-<M)eikrCOs(6'V  .  (2.34) 

o  o 

The  new  path  T  is  thus  the  addition  of  two  Y  paths. 

Similarly,  if  the  above  procedure  is  performed  on  the  second 
boundary  condition,  Eq.  (2.31),  integrated  on  the  path  shown  in 
Fig.  2.6,  with  a  change  of  variable  a'  *  3+a-fT,  the  resulting  equation 
is 


/pP(-0+a+ir) [sin9  -  sina]eikrcosada 


-  /oj.  **  ikrcos(0+(J>  ) 
-2p  (sm-:  -  sin(S-Hj)  ))e  o 

o  o 


(2.35) 


In  addition  to  the  above  integral  equations  there  are  two  more 
conditions  that  are  to  be  satisfied,  namely  continuity  of  pressure 
and  velocity  at  y  ■  0  between  the  two  regions  y  ^  0.  The  two  regions 
are  joined  together  at  <5  ■  0°.  Then  the  continuity  of  pressure  at 
0  ■  0°,  using  Eqs.  (2.14)  and  (2.22)  is  as  follows: 


f  P(cos  i)e  *krcoSada  +  p  e 
Y1  0 


ikrcosO 

o 


_ %  -ikrcosa.  .  _ 

cos  Oe  da  +  P  e 

o 


ikrcos;Q 


(2.36) 


where  the  paths  7^  and  y9  may  be  taken  as  one  of  the  many  possible 


Y-paths . 
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The  transverse  component  of  the  velocity  for  the  incident  wave  is 
given  by 


vi$ 


i  3_ 
kpcr  3<p  p* 


Po 


sin(0-$n)  eikrcos(<j)-<{i0) 
pc 


(2.37) 


and  the  velocity  for  the  diffracted  pressure  is  given  by 

v  »  - — i —  -2 —  p(r,<p)  *  /  P(cosa)sin($+a)e  *krcos($+a)jo  (2.38) 

kpcr  ,r  pc  Jy 

Then  the  condition  for  the  continuity  of  velocity  at  <p*0°  becomes 

/  P(cosa)sinae  ^^rcosotdo  *  -/  P(cosa)sinae  ^rcosa<jot  (2.39) 

T1  Y2 

The  paths  of  integration  for  Eqs.  (2.36)  and  (2.39)  are  given  in 
Fig.  2.7. 

2.5  Solution  for  Plane  Wave  Incidence 

To  solve  the  problem  one  must  attempt  to  construct  a  spectrum 
function  P(coso),  that  will  satisfy  the  boundary  conditions  and  the 
continuity  requirements. 

The  construction  of  the  function  can  be  made  by  closing  the  T-path 
in  Eqs.  (2.35)  and  (2.36)  within  a  region  of  convergence  and  invoking 
the  Cauchy  integral  formula.  If  one  succeeds  in  closing  the  path  T, 
then  the  encircled  poles  would  have  residues  that  satisfy  the  right 
hand  side  of  the  boundary  conditions. 


Fig.  2.8  shows  how  a  T-path  may  be  closed.  The  integrals, 
Eqs.  (2.35)  and  (2.36),  around  a  closed  T-path  must  have  the  sum 
of  their  residues  equal  to  the  right-hand  side  of  these  equations. 
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The  major  characteristics  for  a  function  that  will  satisfy 
Eq.  (2.35)  are: 

1)  P(g-a-TT)  (sin9+  -  sina)  is  even  in  a, 

2)  P(S-a-TT)  (sin6+  -  sina)  has  simple  poles  at  ±(B-<j>o)±2nT 

with  corresponding  residues 

^o  /  .  -+  •  /o  a  w  ikrcos(8-$  ) 

-  (sino  -  sin(S-<£Q))e  ro' 

\ 

Likewise,  the  integral  Eq.  (2.36) 


t/r+/TT  +/..  ]P(-8+-a+7r)  [sin9  -  sina]eikrcosGtda 
1  U1  u2 

-  -2po(sin9'  -  sin(S+^>0))eikrcos(3+4>o)  (2.42) 

suggests  the  following  requirements. 

1)  P(-S+:«+Tr)  (sin9  -  sina)  is  even  in  a, 

2)  P(-3+oc+-~)  (sin9~  -  sina)  has  simple  poles  at  t(B+<}>0)  ±2na 

and  corresponding  residues 


wo  f  .  1 ikrcos(S+0  ) 

-  -rzr  (sin9  -  sin(o+5  ))e  o 

A  1  O 


To  accommodate  all  of  the  above  requirements  the  spectrum 
function  is  considered  as  the  sum  of  two  parts.  One  part  has  poles 
whose  residues  satisfy  the  first  boundary  condition,  the  other 
part  has  poles  whose  residues  satisfy  the  second  boundary  condition. 
Two  such  functions  that  have  the  proper  singularities  and  provide  the 
proper  residues  are  chosen  as 
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If  one  Cakes  Che  firsC  of  Che  above  funccions  and  encloses  the 
simple  pole  at  ,  the  resulting  residue  is 

isini*  -  sin(3-4,o))eikr‘:°s(e-'t0> 

Similarly,  the  residues  at  a»8+4>o  is  given  by 

-  -2-  (sine"  -  sin(3+<}>  ))eikrcos (0+V 

27i  ° 

The  two  functions  happen  to  be  what  is  required  to  satisfy  the 
boundary  conditions.  The  condition  of  evenness  can  be  easily 
verified  by  substituting  3-a-m  and  -B-Kt+m  in  Eqs.  (2.40)  and 
(2.42)  respectively. 

One  further  heuristic  requirement  is  that  for  an  absorbent 
wedge,  solution  of  the  far-field  pressure  must  cancel  out  identically 
at  both  boundaries  due  to  the  absorption  of  the  energy  by  the  wedge 
surfaces.  By  introducing  two  new  functions  that  multiply  Eqs.  (2.43) 
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and  (2.44),  the  final  form  of  P(cosa)  is  obtained. 

P(cosa)  =  M^cOi/^Ca)  -  M2(ct)^2(a)  (2.45) 

where  M^(a)  and  M2(a)  are  the  Maliuzhinets  angle  factors. 


COSTTy .  /28 

M  (a)  *  - r — “ - - -  (2.46) 

sin2^(a+Tr)-sin^o 

cos~o  /28 

M2(  ) - - - 2 - —  (2.47) 

s  in2g(a-Tr ) -s 


and  0^(a)  and  ^(a)  are  defined  by 
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...  .  lo_  [slne~-sin(g-a)l - 

*2  461  [sl„e-+2Sln(^)Sln(^)) 

28+a-H>  -it  — ot-K{>  +tt 

{  sin9++4  tan  ( - — )  tan  ( - ~ — ) 

. ,  2B+a+4>  -it  -a+4>  -ir  _D ,  . 

,  .  ,8+gv  ,  o  ^  /  o  v,  ,  ,B+ok  ,3B+*V  nil  / 

[  s  in  (r~2~)  tan  ( - j - )  tan  ( - £ - )  +s  in  (-£-)  cos  ( — — )  cos  8  ]  >  / 

28+g-H)  -fr  -orHJ>  +ir 

(sinQ  +4tan( - ^ - )tan( - - - ) 

.  .  ,S+cu  26+a+4>0-^  -a+*0+7r  g  B-<P 

[sin(-t~-)tan( - - - )tan( - £ - ) -s in (——) cos — )sin( — ^— )  ]  }  (2.49) 


When  one  of  the  surfaces  is  rigid,  a  special  form  of  P(cosa) 
is  needed  since  there  is  no  absorption  at  that  boundary.  For  these 
special  cases,  the  functions  '^(a)  and  ^ can  be  modified  as 
follows  when  6  «  0° 


Vx(a) 


£o_ 

4Bi 


fsin?*-sin(S-ct)  ] _ 

(sin0++2sin(^^)sin(^z|a)  ] 


S-g  2B-5-0^ 

(sin  (-^-)+2tan( - 


a-0o+7T 

■)  tan  ( - — )  cos 


^  _  2c-.v-:  g— c„+~  23-g- : ,  . 

{sin  O2^) -4  tan  ( - ^ — )  tan( - | — )sin( - £ — ~)cos( — ~)  ; 


(2.50a) 
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,8+cu 
-cos  (— j-) 


^  '  46‘ 


2g+oc+4)o-7r  -oH>0+tt 

(sin9++4tan( - £ - )  tan( - ^ - ) 

Q.  28+a+6  -it  — ct+-4>  +tt  „  / 

[  sin  (-^)  t an  ( - j - )  tan  ( - £ - )+sin  (^)  cos  (-^— )  cos  S  ] }/ 

26+a+6  -7T  -art-<|>0+TT 

{sinQ  +4tan( - r - )tan( - £ - ) 

28+0+6  -TT  — Ot+(J)  +TT  Q.  -Q.  8-6. 

[sin(^y^)tan( - - — )tan( - j - )-sin(— )cos( — —)sin (■-—-) ] }  (2.51a) 


and  when  9  =0° 

,8-cu 
p  -cos  (— j-) 


{ s in9~+4 tan ( 2 - )  tan(a 


„  28-a-6  -it  a-6  +7t  „ 

[sin(^)tan( - )tan(  — f— )+sin(^)cos(^|:::i)cos8] ! 


28-o-60-t  a-60+?r 

(sino  +4tan( - r - )tan( - £ — ) 


2S-a-5.-~  a-6  ++ 


[sin(^j^)tan( - ^ — )tan( - £ — )-sin  (— ^-)cos( — — )sin(— 


)-sin(^)cos(^)sin(-T— )]  (2.50b) 
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,1.  fft-v  .  P°  [sin8~-sin(g+a)1 _ 

2  4ei  [  sinQ  “+2sin(^)  sin  (^i)] 


28+a+6  -77  -a+4>  +rr  / 

{sin(^p)+2tan  ( - - - — )  tan  ( - ^ — -)cos8  )f 


28+a+'t>n-7T  -a+<})  -ht  28+a+<*>  a+6 

(sin(=y^)-4tan( - £ - )tan( - £ - )sin( - - - )cos(— £ — ) 


(2.51b) 


It  is  necessary  to  show  that  the  function  P(cosa),  when  substi¬ 
tuted  in  Eqs.  (2.36)  and  (2.39),  represents  continuous  pressure  and 
velocity  fields.  Eq.  (2.36)  clearly  implies  that  the  pressure  is 
continuous  at  ^  *  0*.  Continuity  of  velocity  requires  further  attention. 
Consider  two  paths  of  integration  in  the  domain  -n  £  a  <_  tt,  where  there 
are  two  possible  poles:  28  -  -  77  and  -tr  +  (!>0*  A.  small  absorption 

should  be  added  to  avoid  wave  generation  at  infinity,  so  that  the 

simple  poles  become  23-6  -  77  -  in  and  -tt  +  <J>  +  in  (Eig.  2.9). 

o  o 

Then  there  are  no  poles  left  in  the  convergent  regions  and  hence  the 
residues  are  zero.  Therefore,  the  continuity  of  the  velocity  is 
identically  satisfied. 

2.6  Far-field  Solution  for  a  Plane  Wave  Incident 

Frequently  in  wave  propagation  problems  the  field  pressure  at 
large  distances  from  the  source  is  desired.  Consider  the  integral 
in  Eq.  (2.14)  with  the  upper  sign  for  the  case  when  0  <  *  < 

,  .  -ikrcos( ‘-ct) 

p  *  i  P(cosa)  e  da  (2.52) 


If  the  non-dimensional  distance  kr  can  be  taken  large  enough 
then  it  is  possible  to  get  an  asymptotic  value  of  the  integral  by 
the  steepest  descent  method.  The  method  is  one  of  finding  the  saddle 
point  (or  extremum  point)  and  integrating  along  the  exponentially 
decaying  steepest  descent  path  which  passes  through  the  saddle  point. 
The  path  of  integration  should  be  chosen  such  that  Re  (-ikrcos  (<f>-ot) ) 
has  a  maximum  value  at  some  point.  Therefore  the  function  has  an 
extremum  if 

(-ikrcos (<})— ot) )  ■  0  .  (2.53) 

Then  the  saddle  points  are  located  at 

a  *  <(>  ±  2nir  (2.54) 

and  the  corresponding  steepest  descent  paths  are  defined  by: 

Im(-ikrcos (0-ci))  *  -kr  *  const  (2.55) 

or 

cos(:-'ir)coshi^  *  1  (2.56) 

The  above  expression  is  the  equation  of  the  steepest  descent 
path  which  passes  through  points  a  *  5  -2n~.  It  is  helpful  to 
know  the  slope  at  the  saddle  points  and  asymptotes  as  ±  00 .  The 
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slope  and  the  asymptotes  are  determined  from  Eq.  (2.56)  as  follows: 


±  coshct.  =  ±1  at  a.  =  0 

l  i 


(2.57) 


cos($-oi  )  =0  as  a±  =  ± 


(2.58) 


ar  =  ^  ±  (2n-l)  y  . 


(2.59) 


For  the  observer  angles  0  £  $  £  8,  the  steepest  descent  path 
is  sketched  in  Fig.  2.10  by  the  use  of  Eqs.  (2.54),  (2.57)  and  (2.59). 

After  completely  defining  the  saddle  point  and  the  steepest 
descent  path,  an  approximation  of  the  integral  in  Eq.  (2.52)  is 
written  as  [31], 


/  P  (cosa)  e-lkrc<JS<;'')d, 
SDP 


/2~  ir./U  -ikr 
'  - —  e  e 


e  P(cos:) 


(2.60) 


It  should  be  noted  that  the  solution  in  Eq.  (2.60)  is  valid 
for  large  values  of  kr  and  an  observer  angle  :  not  close  to  the 
incident  and  reflection  shadow  boundaries  where  the  function  P(cos^) 
becomes  unbounded. 


wS 


immiw 
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2.7  Geometric  Acoustics 

The  acoustic  scattering  by  a  wedge  can  be  presented  in  two  parts, 
a  diffracted  part  and  a  geometric  part.  The  steepest  descent  path 
contributes  the  diffracted  part  which  is  given  in  Eq.  (2.60)  for 
large  kr.  While  the  pole  that  is  surrounded  by  the  closed  path 
contributes  the  geometric  part. 

Consider  the  case  where  the  receiver  angle  <p  is  in  the  region 
26  -  <J>Q  -  tt  <  <j>  <  0.  The  path  can  be  closed  by  the  y-path  and  the 
steepest  descent  path  so  that  the  pole  is  surrounded  as  shown  in 
Fig.  (2.11).  Then  the  solution  in  this  region  becomes: 


,  ,4.  \  Po„  ikrcos(4>-4>  )  ,  , 

p(r,<J)  =  -y  Cr  e  o  +  Pd ( r , v ) 


(2.61) 


where  p^  when  kr  is  large  is  given  by  Eq.  (2.60)  and  is  the 
reflection  coefficient  defined  as: 

sin(S-4>  )-sin@+ 

C  -  - 2 - -  .  (2.62) 

sin(3-0o)+sin5 


If  the  receiver  is  placed  at  0  >  <5  >  23  -  -5  -  ~  (see  Fig. 

2.12)  the  solution  is  given  by: 


/,  .n  _  io  ikrcos(i-0  )  ,  .s 

p(r,,)  =  - —  Cr  e  o  +  pd(r,.) 


(2.63) 
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The  change  in  sign  between  Eqs.  (2.61)  and  (2.63)  represents  the 
sudden  transition  across  the  reflection  shadow  boundary.  Thus,  the 
total  field  in  the  different  regions  shown  in  Fig.  2.13  can  be  written 
as: 

pt  “  pd  +  pi  +  pr  Region  (I) 

=  p.  +  p.  Region  (II) 

d  x 

*  p^  Region  (III)  (2.64) 


CHAPTER  III 


SOLUTIONS  VALID  FOR  MID-RANGE  AND  AT  SHADOW  BOUNDARIES 


3.1  Mid-range  Solution  for  a  Plane  Wave  Source 

In  the  asymptotic  solution  for  the  far-field  there  are 
discontinuities  near  the  shadow  boundaries,  because  the  coefficients 
Mj(a)  and  M2(a)  are  singular  at  these  angles.  These  regions  are 
referred  to  as  transition  regions.  To  have  a  continuous  total  field, 
one  has  to  derive  a  more  exact  solution  which  is  bounded  across  these 
transition  regions.  A  method  was  first  proposed  by  Pauli  [18]  to 
achieve  this  Cask. 

First,  the  spectrum  function  P(cosa)  in  Eq.  (2.43)  can  be 
manipulated  to  have  the  following  convenient  form. 

P(cosa)  »  Mi(a)yi(a)  -  M2(a)'i'2^a) 

«  {[cot  (n  +  (a-<p0))  +  tan  (it+(a+v0))  Jvi(a) 

+  [cot  ^  (ir-(a-$0))  +  tan  (n-(a+p0))  JV2(“)>  (3.1) 

The  integral  for  the  pressure  can  now  be  treated  as  four  separate 
integrals,  each  of  which  contains  one  simple  pole.  The  first  of  these 
can  be  written  as 

Pl  *  j  iy  V'i(a)cot  (it +(a-?0))e  da  (3.2) 

By  replacing  a  by  j-a  and  changing  the  limits  of  integration 
correspondingly  (Fig.  3.1)  the  dependence  of  the  exponential  on  the 
observer  angle  can  be  removed. 
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PI 


1  r  ,  /  -ikrcosaj 

2  JY.'i'l(a"<(>)cot( - Jn - '  6  d 


where 


n 


26 

ir 


(3.3) 


(3.4) 


The  singularity  of  the  integrand  can  be  separated  by  extracting 
the  singularity  explicitly,  i.e. ,  multiplying  and  dividing  the 
integrand  by  a  factor 

cosa  +  cos(2nnN  +  ($“$<>))  (3.5) 

Then  Eq.  (3.3)  becomes 

1  n+($-a-t?0) 

Pi  ■  j  - ) 


cosa  +  cos(2n«N-+(,-a-p0»  _ikt 

- - -  e  da  (3.6) 

cosa  +  cos(2nrrN  +(#-a-<p0)) 


where  N~  is  an  integer  which  nearly  satisfies  the  equation 


2nxN"  -  ($-a-$0)  ■  -  n  .  (3.7) 

This  condition  ensures  that  only  one  pole  is  generated  and  the 
other  poles  of  the  periodic  functions  of  Eq.  (3.5)  are  excluded. 

In  the  method  of  steepest  descent  for  the  far-field  approximation 
mentioned  in  the  previous  chapter,  the  saddle  point  is  taken  as  the 
main  contributor.  The  solution  valid  in  the  mid-range  and  at  the 
shadow  boundaries  consists  of  evaluating  the  integral  along  the  entire 
length  of  the  steepest  descent  path. 
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Letting  the  conformal  transformation 


ikr(cosa-l);  x 


nr. —  -iit/4  .  a 
2kr  e  sin  j 


(3.8) 


so  that 


2.  iff/4,  a 

kr  e  /cos  2 


(3.9) 


and  substituting  Eqs.  (3.8)  and  (3.9)  into  (3.6)  with  limits  to  ®, 


it  follows  that: 


Pl  ■  i/  Tp-  eilT^  e  ikr  /  4»i(({»~a)cot(- 


w-K<p— at— «p0) 


cosa+cos(2nirN”+(^-a-<j)0))  -x2 

- : -  hm  dx  (3.1°) 

x2+ikr(cos(2nnN  +($-a-<j>0)))  Sa 

By  expanding  the  integrand  in  a  Taylor  series  about  the 

transformed  saddle  point  a  »  0  and  retaining  only  the  leading  term,  Eq. 

(3.10)  reduces  to 

■  aa  —X2 

P1  -  i/^|  ein/4  e"ikr  4»i(t>)a"(6-(?o)  /  - S - dx  (3.11) 

"*  x2+ikra  (?-<?0) 


where 


2nu  N~+  Cv>-P0) 

a  (p— P0)  ”  2cos2( - - — ) 

2 


(3.12) 


The  real  integrand  in  Eq.  (3.11)  can  be  expressed  in  terms  of  the 
complex  Fresnel  integrals  (26|. 


+  ib2 


dT  -  2/7  F(b) 


(3.13) 


.  •»  *  (k  *.•  .  .  •  -•  -»  -  -  ■ 
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where 


F(b)  -  eil>2  J  e“it2dT 


Another  form  of  Fresnel  function  Is  given  by: 

b 

f 

o 


F0(b)  -  etb2  /  e'ix2  dt. 


Since 


-it2  .  -  r*  -it2  .  J—  -llt/4 

i  dr  «  2  J  e  dt  »  /tt  e  , 
o 


and 


F(b)  +  F0(b)  -  j  ✓?  e“i1,/4  eib2  , 


then 


F(b)  -  j  /ff  e~ln/4  fiib2  _  Fo(b). 

F(b)  can  be  expressed  in  terns  of  C(b)  and  S(b) 
F(b)  -  ^  etb2  e-1,/4 

(1-/2  .‘"^(CUl-lSlb))] 


where 


C(b) 


/—  /  cosx2dt. 


and 


S(b) 


/  sint2dt. 


Using  Eqs.  (3.19),  (3.20)  and  (3.21),  Eq.  (3.11)  can  be 


as 


(3.14) 

(3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 

rewritten 
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P1-i  I  ♦,<♦)  cot(*-p— )  e‘lkr  eikra"<»-»») 

✓  2  1  2n 

/T^)  (1-/2  eilT^4  (C(/kra”(<{)-90))  -  iS(/kra~(<>-<)>0))) 

(3.22) 

Similarly,  expressions  can  be  written  for  the  remaining  three 
terms  of  Eq.  (3.1).  Then  the  total  diffraction  field  which  behaves 
smoothly  in  the  transition  regions  is  given  by  : 


.  tt  -ikr 
p  -  i  - e 

/7 

■"+($*-<l>o) 

{<i'l(<J>)lcot( - )F1(b-<(>0) 

2n 

TT+((||-hj»o) 

+  tan( - )F*  (<j>+9o“2  S)1 

2n  + 


where 


*-($-<Po>  , 

+  ^(❖Mcot^ — - )f  ‘(r?0> 

+  tan( - ^ - )F*_(4>-H>0+2S) )} 


fJ(x)  -  eikra_(x)  /a_(x) 


(3.23) 


(1-/2  ei7r/4  (C(/kra2:(x))-iS(/kra:(x))) 


(3.24) 


-/  \  n  2/2nwN“-x. 
a  (x)  -  2cos4( — ^ - ) 


(3.25) 


in  which  N4  are  the  integers  which  most  nearly  satisfy  the  following 


conditions  (17) 


50 


2nnN+-x  *  n 


2nirH~-x  *  -7f 


(3.26a) 

(3.26b) 


3.2  Mid-range  Solution  for  a  Line  Source 

The  diffraction  field  generated  by  a  line  source  can  be  written  by 
integrating  the  plane  wave  spectrum  over  two  dimensional  space,  i.e., 


^ 

p(r.*)  -  -S-_  fyfy  P(«,0)e'lkrcos(*^Villr°cos(*°_<’),lc,da  (3.27) 


Using  the  transformation  $  -  a  to  a  and  <f>0  -  a  to  a  and  changing 
the  limits  correspondingly,  one  obtains: 

p(Cl„.£^  If  P(,-a,,0-.)e“llurcoso  (3.28) 

JlT  '  1t2 


Consider  the  second  term  of  the  modified  form  of  P(cosa)  from  Eq. 


(3.1) 


P2  (a, a)  -  j  tan  -jrg  (n  -Ka-o)Hi(a.o) 


(3.29) 


where  $0  is  replaced  by  o  and  n  denotes  2S/*. 

Substituting  Eq.  (3.29)  into  Eq.  (3.28)  leads  to 


e~in/4  7r+$>-a+00-o 

p  .  -  — •  /  /  1>i($-a,$0-a)tan( - j- - )  e“ 

Z  iSS  Ti  t2 


.  -ikrcosa  -ikrQcosa ,  . 

)  e  e  u  doda 


(3.30) 


Then  multiplying  both  the  numerator  and  the  denominator  by  the 


factor 


cos^(  7 +9  o"0“a+ 2ms  N+ )  co s-~(  ( 0 -a+o +2 n^  'i* ) 
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Eq.  (3.30)  becomes 
-iir/4 

„  .  ® _ 

2  4/2? 


*yjy2{  ♦+^0“a“h 

*  *  r  ~  /  \  _ 


[cos( - 2 - )-sinj} 


^2 (|f~a*(i>o  c)  -ikrcosa . -ikr„coso  ,  r_  -n 

- --~~h -  >  e  e  doda  <3*31> 

(cos  ( - 7j - )+sinjj 


where 


h  »  28  +  2nirN* 


(3.32) 


,  rr+tp~a+<p0-o 

<l>l(<i>-a,<t>o-o)  «  ij>1(<{>-a,<jo-0)tan( - ^ - ) 


<p+$0-a~o-h 


(3.33) 


In  the  first  term  of  the  Eq.  (3.31)  changing  a  by  -a,  then 


-in/4  i|>j  ($-h*,$0-o) 


4/2? 


fyfy  t 

Tl  Y2 


(cos(-^^SZil)-sin|j 


9+90-a-h  "~ 

(cos( - 2 - )+sinj] 


,  -ikrcosa  -ikrQcoso .  , 

}  e  e  °  doda 


(3.34) 


■S. 
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where 
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Rl  =  r0  +  r 


(3.41) 


yields 


p2 


12e~i7r/A  e-ik(r+rQ) 
✓2tt 


4>-H|>0— h 

,  <J>0>c°s2  ( - - — ) 


/  pK(p)  e"kRlp2<ip  (3.42) 


with 


2ir  dy 

K(p)«  J  - 

o  _  _  4>+$>o-h  _  $-<?>Q-h 

p2[/r/r0cos2Y+^r^7rsin2Y+sin( - ^ — )sin2Yj+2i/rrQ/K£Cos2( - - — )]  J 

(3.43) 

K(p)  can  be  integrated  by  putting 

z  -  eiY  (3.44) 

so  that 

z  dz 

K(p)  -  -41/  — -  (3.45) 

UNIT  Az“  +  2Bz2  +  C 
CIRCLE 

where 

,  r= — 

A  -  p2{/  f-  -  -  12sin( - = - )J  ,  (3.46) 

rn  r  i 
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Again  introducing  a  new  variable  A 
iA2  »  ik(Rj-S)  +  kRjp2 


(3.52) 
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Then  the  integral  in  Eq.  (3.50)  becomes 


«  ~ix2  ik(RrS) 

*  iX/kRi  e  e  1 

/  - : - - -  dX 

_  v/iX2wiX2  S__  , 

/icTRpsT  * 'kRj  kRj  2i  Ri  } 


f  eik(Rl_S)  -iX2  ^ 

j  ________  e  dX 

/k(B."i“S)  /X2  +  2Sk 


(3.53) 


A  Fresnel  approximation  of  Eq.  (3.53)  can  be  obtained  when  kRj  »  1 
by  replacing  the  non -exponential  part  of  the  integrand  by  its  value  at 
the  lower  limit,  to  give  the  approximation 


eik(Rj-S)  j  e~ix2dX  - - i -  F(/k(R1-S))  (3.54) 

✓kUi+S)  ✓k(R1-S)  /k(Ri+S) 

where  F(b)  is  given  in  Eq.  (3.19). 

Finally,  the  integration  of  Eq.  (3.30)  is  reached 


4ie~iir/VikRlj>l(»t(?Q)—  1  -  F(/k(R!-S)) 


/k(Ri+S) 


2  —  V ,  (  5  *  <J>  o  )~ 
/2 


-ikS 


(1-/2  eilT/4(C(/k(R1-S))-iS(/k(R1,-S))) 


«^k(Rj+S) 


(3.55) 


with 

and 


h  -  23  +  2nnN'1 


/  9-^o-h 

S  •  /(r+rQ)2  -  4rrQcos2( - ~ — ) 


2 
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Similar  expressions  can  be  found  for  the  ocher  three  terms  of  the 
Eq.  (3.27).  Then  the  total  solution  for  a  line  source  when  k(r+rQ)  »  1 
is 


7T+(<p-^0)  2  *+(<f+lpo)  2 

P  m  n{'?j(6.'}>o)Icoc( - — )  F_(<J>-<f>0)  +  tan( - ^ —  )F+  ( <f>-Np0— 2B )  ] 

1T-($-4>0)  2  w-(<p-Hp0> 

+  ^  o  H  cot  ( - 2n - ^F+^',<i>o)  +  cat»(  - 2n - ^F-  ^'hS’o+2B^^ 

(3.56) 

where  2  -ikS?(x)  r-r-j~ 

F-;(x)  *  -e  - — -  /a~(XJ 

✓kCRj+S^x) 


(1-/2  e^^tf/kCRj+SgU))  -  iS(/k(R1+S*(x))  ) 

(3.57) 


S*  *  /(r+r0)2  -  2rr0a± 


(x) 


(3.58) 


and 


.=(*>  -  2coS2 


(3.59) 


is  given  in  Eqs.  (3.26a)  and  (3.26b). 


3.3  Mid-range  Solution  for  a  Point  Source 

The  diffracted  pressure  for  a  point  source  can  be  written  in  the 

form 


P 


'  Y 


P(cosa) 

kRUl 


e~ikR(a)da 


( 3.b0) 


where 


R(a)  -  /r2  ■+  +  2rr0  cos(y-a>  +  (z-zQj2 


(3.61) 
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and  P(cosa)  is  defined  in  Eq.  (2.45). 

The  function  is  not  single  valued,  so  the  branch  cuts  are  located 


R(o)  =  0 


(3.62) 


i  t  <2  .  , 


r2  +  r2  +  (z-zQ)2 


-  cos(ij)-a) 


-  cos(({)-ar)coshOi 


(3.63) 


Since  (r-rQ)2  =  r2  +  r2  -  2rra  >  0,  then  r2  +  r2  >  2rr0  which 


implies  that 


r2  +  r2  +  (z-z0)2| 


>  1  . 


(3.64) 


As  a  result,  the  branch  points  are  located  by 


cos($-ar)  <  1 


(3.65) 


coshct£ 


r2+r2+(z-z0)2 


(3.66) 


ar  -  (2n-l)ir+<> 


(3.67) 


*  -  cosh~U- 


r2+r2+(z-z0)2 


(3.68) 


The  regions  of  convergence  of  the  integral  in  Eq.  (3.60)  are 


determined  by  the  real  part  of  the  exponential  power. 


58 


Re(-ikR(a))  <  0  (3.69) 

or 

Im(R(a) )  <  0  (3.70) 

The  function  R(a)  can  be  written  in  polar  form  as 

R(a)  »  |R| 1/2  e*^  *  (cosij;  +  isiny)  (3.71) 

From  Eq.  (3.70)  the  condition  of  convergence  reduces  to  siniji  <  0. 
Then  the  regions  of  convergence  are  governed  by: 

it  +  <  ctr  <  <j>  when  a*  <  0,  (3.72) 

and 

it  +  $  >  QLr  >  $  when  >  0.  (3.73) 

From  Eqs.  (3.67),  (3.68),  (3.72)  and  (3.73)  the  branch  cuts  and 

regions  of  convergence  are  illustrated  in  Fig.  3.2. 

To  obtain  the  solution  to  the  integral  in  Eq.  (3.60),  a  procedure 
is  employed  similar  to  the  plane-wave  and  line  source  incidences.  The 

solution  is  divided  into  subintegrals  and  poles  are  extracted 

explicitly.  First  of  these  subintegrals  (using  Eq.  (3.1))  is 
j  TT+(a-90)  -ikR(a) 

PI  -  2  ‘i  *l(a)  cot< - 2^—- >  t  KU)'  do-  <3-74> 

One  can  avoid  the  singularity  again  by  multiplying  and  dividing  the 
integrand  by  a  factor 

2nirN  -  (o-?0) 

cos  ( - 2 - )  (3.75) 

Consider  closing  the  y-path  along  ;  +  i»  to  ;  -i<*.  Then  p^  is 
given  by 


J 


Fig.  3.2  Regions  of  convergence  and  branch  cuts  for 
point  source. 


'.1 


m 


00 

P1  “  "  J  ^  ((4>+i<*i) 


2nirN"  -  (<p+iai~90) 
cos  ( - = - ) 


iH3c/r2+ro+2rr0coshai+(z-zo)2 


k/  r^+r^+^rr Qcosha j+(  z-zQ  j 


(3.76) 


m 


where 


<J>j  (<J>+iai)  - 


n+6+ia^-(}i0  2nitN  -($+ia^-<>0) 

<|>j(  <►«<*!  )cot( - ^ - )  cos  ( - ^ - ) 


Using  the  fact  that  pj(a)  ■  pi(-a) 
00 

Pi  -  -i  /  1 ( $+iai  ($-!<*£  )] 


(3.77) 


,(2---  1--°-)  cosh  %  e-ik/r2+ro+2rr0coshai+(z-Zo)2 


[coshaj+cos(2nirN  -(4>-$0))J  / r^+r£+2rr0coshui+(  z-zQ)‘ 


(3.78) 


An  asymptotic  expression  for  the  integral  in  Eq.  (3.77)  can  be 
obtained  by  letting 


^ rro  -i.it/ U 
t  *  — ^ -  e  sinhai/2 


(3.79) 


where  Rj 


+  (z-z. 


R(a)  -  /(r+r0)J  T  U-3,^  (!+ 


-  Ri  /1-iT- 


nmwstuwWi  il*  iviv  WKFsynr*  •.  *  r  -.-■  v.  *.-  •--  *--•  --■ 


TT: 
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Rewriting  Eq.  (3.78)  with  new  variable 


m 


pi 


b  _iir/4  •  fiP1(4>-*-ioii)+^1(<p-iQti) J  2n7rN 


kRf 


/ 


cos(- 


✓  1  -  it2 


^-ikRi/I-ir7 


4rr0  2nirN“-(<ji-<t)Q) 


i2+i  — = —  cos2( 


-) 


dr 


(3.81) 


Expanding  t'}'i(4'+i<*i)+4'i($-ioii) J//  1-ix2  in  a  Taylor  series  about 
0  and  retaining  the  leading  term,  Eq.  (3.81)  becomes: 


PI 


Uftr —  *  2n1tN  -<*-*o) 

_ °  e“in/4  ’J'jU)  cos( - j - ) 


kR, 


-ikRj/l-it2 


dt 


4rru 
x2+i — —  cos2( 

*1 


o  2nirN  -($-<|>0) 


2n 


(3.82) 


In  order  to  integrate  the  integral  in  Eq.  (3.82)  consider  the  case 
kRj  >  1,  then  /l-it2  =  l-it2/2  and 


„  „  2  . -in/4  -ikRi  J2rr°  ” 

Pi  *  '^(ijOe  e  1  / — z~  a  (*) 


a-kRLT2/2 


r2  +  i  ~VTr 

\ 


dr 


t4  +  i  -=^a  a  (x) 

R. 


(3.83) 


with  x  ■  q>— <p0 
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where 


.-(X)  -  2co,2<i2I»-2‘) 


(3.84) 


which  has  a  similar  form  to  Eq.  (3.11). 

Defining  Fresnel  integral  of  the  form 

o»  -tT2 

bj  £ -  dt  -  2/?  F(b/tf) 

-  x2+ib2 

and  using  Eqs.  (3.13)  through  (3.15)  the  final  form  is  obtained: 

,  2ir  , *+($-$0)  -ikRi  ikrrn/Ri  * 

Pi  *  1  — ~ -  <K$)cot( - ^ - )  e  1  e  °  1  a  (x) 


(3.85) 


/2kRj 


2n 


/a~ (x)  (1-/2  e1,r/4(C(/ (krr0/R|)a  (x))-  iS(/(krr0/R|)a  (x))) 

(3.86) 


There  are  similar  expressions  for  the  other  three  terms;  then  the 
total  solution  becomes: 


PI  "  «  1  {V^lcot*-^ - >  F 


"+($*»  0>  , 

+  tan( - - )  F^^+^o"20^ 

*-($-$o) 

+  ’^2(<?){cot( - ^ - )  F^($-<?o) 

w-(<?+90) 

+  tan( - = - )  F3(^0+2S)J}  (3.87) 

zn  - 


where 


F*U) 


ikrr0/Ri  ai(x)  j  z.  \ 

•  e  1  / a  (x) 

(1-  /2  ei,T/4(C(/(krr0/R1)ai(x))-iS(/(krr0/R1a-(x))) 


CHAPTER  IV 


NUMERICAL  RESULTS 

The  solution  for  the  acoustic  diffraction  from  impedance  covered 
wedges  is  a  function  of  several  physical  and  geometric  parameters;  the 
surface  impedances  9*  and  9~,  the  source  angle  $0,  the  non-dimensional 
wave  numbers  kr  and  kr0  and  the  wedge  angle  9.  In  this  chapter, 
diffracted  and  backscattered  pressures  are  plotted  for  different  values 
of  these  parameters.  In  order  to  assess  the  influence  of  one  of  the 
parameters,  the  others  are  kept  constant.  Most  of  the  graphs  shown  in 
this  section  are  polar  plots  of  the  angular  distribution  of  the 
diffracted  pressure  in  decibels  relative  to  the  amplitude  of  the 
Incident  pressure. 

The  mid-range  solution  of  diffraction  of  a  plane  wave  in  Eq. 

(3.23)  diminishes  at  a  rate  proportional  to  /kr  because  of  the  factor 
/kr  in  the  argument  of  the  Fresnel  integrals.  To  study  the  influences 
of  the  wedge  angle  and  the  surface  impedances  on  the  angular  variation 
of  the  pressure,  it  is  advantageous  to  remove  the  geometric  spreading 
factor  (/kr)  from  the  expression  for  the  diffracted  pressure  for  plane 
wave  Incidence.  Thus,  henceforth  only  the  product  of  /kr  •  Pd(r,p)  is 
plotted.  By  the  same  reasoning,  the  diffracted  pressure  for  a  point 
source  is  normalized  by  the  incident  pressure  at  the  edge  due  to  a 
point  source  e“*kro/kr0  an(j  multiplied  by  /kr  in  oraer  to  facilitate 
comparisons  between  the  different  plots.  For  all  the  plots  for  a  point 
source,  the  source  and  the  receiver  are  placed  in  the  z*0  plane. 

Because  of  limitations  of  the  plotter  the  following  symbols  are  used 
for  labeling: 
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kr 

-  KR 

0+  a 

BAT 

$0 

=*  SA 

krQ 

-  KRO 

6"  “ 

BAB 

3 

=  WA 

4.1  Dependence  on  Surface  Impedances 

In  Che  derivation  of  the  solution,  the  Brewster  angles  and  6 
were  used  instead  of  complex  impedance  Z~ .  Rewriting  the  relation 
between  Z±  and  3*  as 


sin0~  =* 


pc 

Zr 


various  impedance  conditions  ranging  from  pressure  release  to  rigid  are 
considered.  The  list  of  impedances  and  their  corresponding  Brewster 
angles  are  listed  below 


Z  (Impedance)  3  (Brewster  Angle) 

Rigid  «  0° 


2.00  pc 

30° 

1.16  pc 

O 

O 

vO 

Matched 

PC 

90° 

Pressure  Release 

0 

90°-i5 

The  influence  of  the  impedance  of  the  insonified 
and  uninsonified  surfaces  are  considered  in  the  first 
These  figures  are  plotted  by  employing  Eq.  (3.23)  for 
parameters: 

1.  Fig.  4.1  presents  results  for  the  diffracted 


(source-facing) 
five  plots, 
the  following 

pressure  for  the 


impedances  *  30° ,  60°  and  90°  on  the  insonified  surface 
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while  6“  ■  0°  (rigid)  remains  the  same  on  the  uninsonified 
surface  with  kr  *  20,  <p0  ■  60°  and  3  *  135°. 

2.  Fig*  4.2  presents  plots  of  the  diffracted  pressure  for  the 
source  angle  $0  *  -60°  to  show  the  influence  of  the  impedances 
of  the  uninsonified  surface.  All  other  parameters  are  the 
same  as  in  Fig.  4.1. 

3.  The  influence  of  an  imaginary  component  of  impedance  is 
considered  in  Figs.  4.3  and  4.4  where  6+  =  90°  -i5  and  9“  * 
0°.  The  angle  of  incidence  is  $0  =  105°  in  Fig.  4.3,  while 
<j0  *  -105°  in  Fig.  4.4  and  kr  *  20  in  both  figures. 

4.  The  diffraction  by  a  wedge  covered  by  absorbent  surfaces  with 
0+  *  60°  and  0~  «  30°  are  shown  in  Fig.  4.5.  The  values  for 
kr  is  20,  $0  is  105°  and  $  is  135°. 


4.2  Dependence  on  the  Wedge  Angle 

To  facilitate  comparisons  of  diffraction  from  different  wedge 
geometries,  plots  for  four  wedge  angles  are  made  by  using  £q.  (3.23). 
Since  the  diffracted  field  is  mainly  controlled  by  the  shadow 
boundaries  and  nulls,  for  each  wedge  angle  two  different  source  angles 


are  used  to  demonstrate  their  influence.  In  Figs.  4.6  through  4.9  the 


angle  of  incidence  is  <?0  ■  110s,  whereas  in  Figs.  4.1U  through  4.13, 

$0  *  30°.  The  values  of  other  parameters  are  kr  ■  20,  6+  *  30°  and  0~ 

-  0°. 


4.2  Diffraction  of  plane  wave  with  impedance  cover 
on  the  uninsonified  surface  with  kr  =  20,  9“  = 
-  -60°  and  3  =  135°. 


JKR-  20  BAT »  90  -51  BAB-  0  SA-  105  WA-  135 


Fig.  4.3  Diffraction  of  plane  wave  with  total  absorption 
on  the  insonified  surface  with  kr  =  20,  =  90 

5“  =  0°,  5  =  105  and  £  =  -35°. 


Fig.  4.6  Diffraction  of  plane  wave  fron  a  120°  wedge  with 
$  *  110°.  0+  =  30°.  6“  =  08  and  kr  =  20. 
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4.3  Dependence  on  krn 

When  che  incident  pressure  field  is  due  Co  a  point  source,  another 
parameter  krQ  is  introduced  which  is  the  direct  distance  from  the  tip 
of  the  wedge  to  the  location  of  the  source.  As  the  source  is  moved 
away  from  the  wedge,  the  incoming  spherical  waves  impinging  on  the 
wedge  look  like  plane  waves.  To  demonstrate  the  influence  of  krQ  on 
the  diffracted  field,  five  plots  are  made  using  Eq.  (3.87)  with 
impedances  0+  *  60°,  0“  =  30°,  kr  =  20,  4>0  =  60°  and  B  ■  135°.  The 
non-dimensional  parameter  krQ  takes  the  values  5,10,  20,  200,  2000  in 
Figs.  4.14  through  4.18. 


4.4  Excess  Attenuation 

It  is  practical  to  model  Che  diffracted  field  from  a  point  source 
as  that  due  to  a  plane  wave  with  a  correction  factor  for  spherical 
spreading.  This  model  assumes  that  there  is  a  geometric-optic 
relationship  between  spherical  and  plane  waves.  Fig.  4.19  is  plotted 
to  explore  this  relationship.  The  veritcal  axis  of  Fig.  4.19  is  in 
terms  of  excess  attenuation,  defined  as 


Excess  Attenuation  ■  -20  logjg 


Pd 

Pdir 


dB, 


(4.1) 


where  p^  is  the  diffracted  field  from  Eq.  (3.87)  and  pd^r  is  the  direct 
field  (e^r2/kf2)  at  the  receiver  in  the  absence  of  a  wedge.  The 
horizontal  axis  is  the  logjg  of  non-dimensional  distance  kr  from  the 
tip  of  the  wedge  to  a  receiver  located  at  *0  *  -9U°  from  the  x-axis. 


Fig.  4.14  Diffraction  of  point  source  radiation  with  kr 
0  -  60°.  9'  -  30°.  i  -  60°  and  2  *  135°. 


Fig.  4.17  Diffraction  of  point  source  radiation  with 
kr  =  200,  S+  =  60°,  =  30°,  :0  =  60°  and 
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4.5  Backscattering  from  a  Wedge 

For  many  applications,  it  is  imperative  to  calculate  the 
diffracted  pressure  at  the  source  location,  the  so  called  monostatic 
diffraction.  The  next  six  plots  exhibit  this  backscattered  pressure. 
For  this  situation,  the  source  and  the  observer  are  at  the  same 
location,  i.e.,  <J>  ■  $0  in  Eq.  (3.23).  In  Figs.  4.20  through  4.23,  the 
backscattered  pressure  is  calculated  for  wedges  whose  surface 
impedances  are  0+  *  60°  and  0-  ■  30°  with  kr  *  20.  In  the  next  two 
plots  the  impedances  on  the  two  surfaces  are  equal  9+  *  0~  =  30°,  while 
kr  -  20  in  Fig.  4.24  and  kr  -  200  in  Fig.  4.25. 

4.6  Comparison  with  the  Haliuzhinets  Solution 

To  conclude  this  section,  plots  comparing  results  of  Maliuzhinets 
and  results  from  Chapter  II  are  presented  in  Figs.  4.26  and  4.27.  The 
plots  are  for  the  far-field  diffracted  pressure  of  plane  waves  derived 
from  Eq.  (3.23).  Two  wedge  angles. namely ,  0  ■  135°  and  0  *  165°,  are 
used  in  the  plots  because,  for  these  special  wedge  angles,  the 
complicated  Haliuzhinets  functions  reduce  to  a  finite  product  form. 

The  resulting  solution  is  adapted  from  Skudrzyk  [7].  In  both  figures 
0+  -  60°,  0”  -  30°,  kr  -  20  and  $0  -  60°. 


KR-  20  BAT-  60  BAB-  30  SA-  110  WA-  165 


Fig.  4.22 


Backscattering  of  plane  wave  from  a  165°  wedge 
with  9  ■  60®,  9  m  30®,  kr  *  20  and  ;  *  110®. 


CHAPTER  V 


DISCUSSION  OF  NUMERICAL  RESULTS  AND  CONCLUSIONS 

5.1  Introduction 

The  primary  objective  of  this  study  is  to  develop  a  solution  for 
diffraction  of  acoustic  waves  by  a  locally  reacting  impedance  covered 
wedge.  This  objective  is  achieved  by  constructing  a  spectrum  function 
from  the  physical  interpretation  of  the  singularities.  The  resulting 
solution  from  this  analysis  is  a  closed  form  solution,  which  is  a 
significant  improvement  on  the  Maliuzhinets  solution. 

In  this  chapter  discussions  of  the  physical  interpretation  of  the 
diffraccion  problem  solved  in  Chapter  II  and  the  numerical  results  that 
were  compiled  in  Chapter  IV,  using  the  mid-range  solutions  from  Chapter 
III,  are  presented. 

5.2  Physical  Interpretation  of  P(cosg) 

The  spectrum  function  P(cosa)  (see  Eq.  (2.45))  consists  of  two 
parts;  the  two  angle  factors  Mj(a)  and  M2(a)  multiplied  by  two 
diffraction  factors  ^i(a)  and  ^(a),  respectively.  The  geometrical 
shadow  boundaries  are  represented  in  P(cosa)  and  are  contained  in  Mj(a) 
and  M2(a).  Although  they  are  different  from  the  Sommerfeld  angle 
factors,  they  possess  identical  poles.  Four  of  these  poles  represent 
Che  following  geometrical  boundaries: 

a.  23  -  >o  ~  T  represents  the  angle  or  specular  reflection  from 
the  upper  surface  (+3). 

represents  the  angle  of  the  incident  shadow 
boundary. 


b.  +$0  -  t» 


c.  -23  -  $0  +  n  represents  the  angle  of  specular  reflection 

from  the  lover  surface  (-3). 

d.  +$0  +  u  represents  the  angle  of  the  incident  shadow 

boundary. 

The  functions  vi(a)  and  ^>2(a)  are  relatively  insensitive  to  the 
argument  a  except  when  a  is  close  to  the  zeros  and  poles  of  P(cosa). 
For  example ,  when  a  approaches  26  -  $0  -  x  ,  representing  the 

reflection  angle  from  the  upper  surface,  approaches  the 

reflection  coefficient  Cr  (see  Eq.  2.62).  Although  ^(a)  and  i|»2(a)  do 
not  have  any  poles  in  the  region  Re|a|  <  6,  they  have  two  zeros  each 
when  the  exit  angle  is  equal  to  one  of  two  Brewster  angles  of  an 
Impedance  covered  surface.  The  nulls  az  of  the  functions  can  be  given 
in  terms  of  the  grazing  Brewster  angles  0+,  and  6~  as  follows: 

For  Vl(<0 

az  “  3  -  0+  (5.1) 

and 


and  for  ^2^a^ 


az  ■  3  -  x  +  0+ 


az  »  -6  +  9~ 


(5.2) 

(5.3) 


and 


az  *  ~3  +  x  -  9” 


(5.4) 


The  function  P(cosa)  has  another  property;  the  term  Hj(u)  yj(a) 
interchanges  with  M2(n)  when  the  signs  of  a,  representing  the 

observer  angle  ?,  and  are  interchanged  as  well  as  j+  by  This  is 

expected,  since  the  chosen  coordinate  system  is  symmetric  about  the 


axis  x-0 
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f sin0+-sin(0-a) ] 


sin6+'+sin(0-a)  1 
P(  20  a;  [sine-+sin(B+a)J 


It  is  apparent  from  Eqs.  (5.7)  and  (5.8)  that  only  wedges  with 
rigid-rigid  or  soft-soft  boundary  conditions  have  spectrum  functions 
P(cosa)  with  periodicity  40. 


5.3  Discussion  of  Numerical  Results 


5.3.1  Transition  Regions 

The  mid-range  solutions  (Eqs.  (3.23),  (3.56)  and  (3.87))  for 
diffraction  of  an  acoustic  plane  wave  consist  of  two  terms;  each  term 
having  two  F+(x)  terms  defined  by  Eqs.  (3.24),  (3.57)  and  (3.88).  At 
the  incident  and  reflection  shadow  boundaries,  the  far-field  solution 
(Eq.  2.58)  becomes  unbounded.  This  discontinuity  in  the  geometric- 
optics  field  is  compensated  separately  by  one  of  the  tour  F*(x)  terms 
in  the  mid-range  solution. 

The  F+(x)  functions  contain  an  argument 


a  (x)  *  2cosn - - - ) 


which  is  a  measure  of  the  angular  separation  between  the  field  point 
and  the  incident  or  the  reflection  shadow  boundaries  which  are  listed 
in  section  5.2.  At  one  such  boundary,  one  of  the  cotangent  or  tangent 
functions  in  the  mid-range  solution  becomes  singular.  The  integers  N* 
and  N~ ,  which  are  defined  by  Eqs.  (3.26a)  and  (3.26b),  assure  that  the 
argument  of  a±(x)  becomes  zero  at  that  boundary.  Since  the  a*(x) 
function  appears  as  a  multiplier  in  F+(x),  the  corresponding  term  of 
the  solution  vanishes  but  the  other  three  remaining  terms  are  bounded. 
Because  the  F+(x)  functions  cause  the  solution  to  become  smooth  across 
these  boundaries,  they  are  called  transition  functions. 

For  example  consider  Fig.  4.1  with  kr  *  20,  0+  -  60°,  0~  *  0°,  <j>0 
*  60°  and  3  *  135°.  The  complex  values  of  pressure  on  either  side  of 
the  reflection  shadow  boundary  at  26  -  $0  ~  ^  *  30°  are: 

pd  -  0.0092  -  i0.0295,  $  -  30°', 

-  -0.0128  +  10.0202,  <p  -  30o+. 

Using  Eqs.  (2.61)  and  (2.63)  the  reflected  pressures  at  either 
side  of  the  reflection  shadow  boundary  are: 

pr  -  -0.0111  +  10.0249 ,  $  -  30°~, 

-  0.0111  -  i0.0249,  $  -  30°+ 

Addition  of  reflected  and  diffracted  fields  produce  a  continuous 
field  as  follows: 

pt  -  -0.0019  -  i0.0046  +  pif  $  «  30°", 

-  -0.0017  -  10.0047  *  Pl,  $  -  30o+. 

If  similar  calculations  are  made  for  the  incident  shadow  boundary 
at  -n  +  ?0  ■  -120°  the  following  continuous  pressure  field  results 
across  that  boundary: 
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Pt  =»  0.1816  -  10.5725,  -  -120°-, 

-  0.1874  -  10.5735,  <j>  -  -120o+. 

5.3.2  Influence  of  Impedance  Coverage 

In  Figs.  4.1  and  4.2  the  source  is  placed  above  and  below  the 
wedge,  respectively.  When  the  impedance  covered  surface  is  insonified 
the  magnitude  of  the  impedance  cover  dramatically  influences  the 
diffraction  field.  As  the  impedance  of  the  insonified  surface  is 
increased  thereby  increasing  the  absorption,  the  diffracted  pressure  in 
the  illuminated  zone  decreases.  However,  the  impedance  cover  of  the 
surface  in  the  shadow  zone  has  very  little  effect  on  the  diffracted 
pressure  in  the  entire  field. 

If  the  upper  surface  of  the  wedge  has  a  pressure  release  cover  and 
the  bottom  surface  is  rigid,  and  the  angle  of  incidence  is  on  the 
pressure  release  side  and  close  to  the  grazing  angle.  Fig.  4.3,  the 
diffracted  pressure  in  the  insonified  zone  is  negligible  because  the 
diffracted  energy  is  primarily  absorbed  by  the  pressure  release 
surface.  The  opposite  is  true  when  the  rigid  surface  is  insonified 
(Fig.  4.4),  the  diffracted  pressure  is  significant  over  the  entire 
space.  Fig.  4.5  demonstrates  how  the  impedance  cover  on  both  surfaces 
influences  the  diffracted  pressure.  Since  the  upper  surface  is  oore 
absorbent  than  the  lower  surface,  the  diffracted  pressure  is  higher  in 
the  lower  region. 

5.3.3  Influence  of  the  Wedge  Angle 

The  diffracted  field  is  controlled  by  the  position  of  the  observer 
with  respect  to  the  incident  and  reflected  shadow  boundaries  and  the 


nulls.  Their  effects  are  demonstrated  in  Figs.  4.6  through  4.13  by 
changing  the  wedge  angle  for  different  source  locations.  For  any  one 
geometric  configuration,  there  can  be  only  two  geometric-optics 
boundaries;  one  incident  and  one  reflected  shadow  boundary  or  two 
reflected  shadow  boundaries.  Figs.  4.10  and  4.11  are  examples  of  two 
reflection  shadow  boundaries.  The  relative  position  of  the  zeros  of 
^l(a)  and  with  respect  to  the  angle  of  incidence  causes  sharp 

minima  shown  in  these  figures. 

5.3.4  Influence  of  the  Sphericity  of  the  Point  Source  Field 

The  diffraction  solution  for  a  point  source  is  given  by  Eq. 

(3.87).  The  parameter  R}  which  appears  in  the  solution  is  the  shortest 
path  the  wave  travels  from  the  source  to  the  edge  of  the  wedge  and  from 
there  to  the  receiver.  It  can  be  seen  that  the  reciprocity  theorem  is 
satisfied  in  this  solution.  If  the  values  for  kr  and  krQ  and  {>  and  $>0 
are  interchanged  in  Eq.  (3.87),  which  is  essentially  interchanging  the 
positions  of  the  source  and  the  receiver,  there  is  no  change  in  the 
solution.  Figs.  4.14  through  4.18  illustrates  that,  as  the  source  is 
placed  farther  away  from  the  wedge,  the  normalized  diffracted  pressure 
field  approaches  to  a  plane  wave  incidence  except  near  the  shadow 
boundaries. 

The  next  plot  (Fig.  4.19)  shows  the  excess  attenuation  (see  Eq. 
4.1)  versus  non-dimensional  distance  from  the  wedge.  *  The  curve  almost 
becomes  a  straight  line  for  observer  distance  kr  over  50,  i.e.,  the 
diffracted  pressure  approaches  its  asymptotic  approximation  of  l//kr. 
This  means  that  the  asymptotic  solution  can  be  used  for  observers 
located  at  least  8  wavelengths  away. 
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5.3.5  Backseat tered  Field 

The  plots  of  the  backscattered  pressure  field  show  clearly  defined 
maxima.  The  maxima  are  located  where  the  reflection  coefficient  Cr 
reaches  its  maximum  such  that 


l~sin9* 
r  *  l+sinS^ 

and 

„  l~sin9“ 
*  1+sine- 


i.e.,  at  4> 


i.e.,  at  $ 


e  -  ir/ 2  , 

+  ir/ 2. 


Also  at  these  locations  the  angle  factors  Mj(a)  and  M2(a)  become 
infinite,  so  the  impedance  cover  only  controls  the  magnitude  of  the 
maxima,  but  not  the  location. 

In  the  backscattering  from  a  rigid-rigid  wedge  there  is  no  return 
if  the  source  is  located  at  $  -  0° ,  whereas  Figs.  4.20  through  4.25 
illustrate  that  the  nulls  are  dependent  on  the  impedance  cover. 


5.4  Conclusions 

In  the  discussion  of  the  closed  form  solution  and  the  numerical 
results  for  the  diffracted  pressure  from  impedance  covered  wedges,  many 
significant  factors  affecting  the  diffraction  are  pointed  out. 

It  is  evident  from  the  plots  that  the  impedance  of  the  surface 
lying  in  the  shadow  zone  has  little  effect  on  the  diffracted  pressure 
while  the  impedance  cover  on  the  insonifiea  surface  has  a  significant 
effect  on  the  diffracted  pressure  in  the  entire  space.  The  impedance 
coverage  of  the  insonified  surface  has  similarly  a  dominant  influence 
on  the  behavior  of  the  backscattered  pressure. 
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The  location  of  the  maxima  and  the  minima  of  the  diffracted  field 
may  be  adjusted  by  changing  the  impedance  of  the  insonified  surface  and 
the  wedge  angle.  It  was  also  shown  that  the  plane  wave  asymptotic 
solution  for  an  observer  in  the  far-field  with  the  correction  factor 
for  a  point  source  is  a  good  approximation  for  kr  >  50  except  near  the 
shadow  boundaries. 

In  the  derivation  of  the  mid-range  solutions,  only  the  first  term 
of  the  Taylor  series  expansion  is  included  in  the  solution.  It  has 
been  previosly  established  by  the  author  that  in  the  integration  for 
the  asymptotic  series ,  the  odd  terms  vanish  and  the  high  order  even 
terms  have  very  little  contributions  except  for  receivers  located  very 
close  to  the  absorbent  wedge  surfaces.  The  first  term  in  the  series 
was  found  to  provide  sufficient  accuracy. 

5.5  Suggestions  for  Further  Research 

In  the  mid-range  solutions,  although  the  first  term  in  the 
asymptotic  series  provides  sufficient  accuracy,  an  investigation  of 
higher  terms  near  the  wedge  surfaces  can  be  studied.  The  higher  order 
•  terms  in  the  Taylor  expansion  contributes  the  surface  waves. 

The  problem  of  diffraction  of  acoustic  waves  from  a  wedge  with 
locally  reacting  boundary  conditions  can  be  expended  to  diffraction 
from  elastic  wedges  made  of  plates  or  solid  bodies. 

Furthermore,  the  solutions  obtained  in  this  study  may  be  adapted 
for  problems  of  diffraction  of  sound  over  a  wedge  located  on  a  rigid 
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APPENDIX  A 


CONDITION  OF  ZERO  VALUE  OF  THE  LOOP  INTEGRAL 

A  loop  Integral  is  zero  [7], 

I  S(z)  e-ik"°s<2>d2  .  0, 

°1«2 

if  the  function  S(z)  is  even,  i.e.,  if  S(z)  =  S(-z) 
where 

z  ■  u  +  iv 

Proof : 

The  path  of  integration  is  deformed  into  lines  arbitrarily 
close  to  u  *  -  T  and  u  *  +  tt  .  Since  cos(z)  =  cos(-z)  and  S(z)  * 
S(-z)  the  integrands  along  the  line  0  >_  v  >  -  »  ,  u  ■  -  ~  and 
0<^v<°°  ,  u  ■  "  are  equal.  But  since  the  directions  of  the  two 

paths  are  opposite  their  total  contribution  is  zero. 
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COMPUTER  PROGRAM  LISTINGS 
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